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' Abstract. In this article we review the quantization of the Dirac-field on a 

' curved spacetime. For that purpose we describe the construction of the local 

observable algebras in the algebraic approach to quantum field theory. Among 
the possible states we single out the so called Hadamard-states, which are the 
ones relevant for physics. Finally, as an example, we give a definition for an 
adiabatic vacuum state of the Dirac-field on a Robertson- Walker spacetime. 
We believe that these states are physical in the sense that they have the sin- 
£\J i gularity structure of Hadamard-form, although we cannot give a formal proof 

of this conjecture. 
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qq ■ 1. Introduction 

| In this paper we are concerned with quantum field theory on a curved spacetime 

. (QFT on CST). Among the vast literature in this area on the Klein- Gordon-field 

00 ' only little work is done for the Dirac-field. The reason for this are not concep- 

0^ ■ tual problems but the greater technical complexity involved in the description of 

multicomponent fields. 

After Hawkings discovery that a black hole not only absorbs the energy incident 
O i' on it, but also loses energy through the emission of pairs of elementary particles 

there were a great interest in the treatment of quantized fields in a classical back- 
ground. This semiclassical description of quantum fields propagating in a non flat 
background can serve as an approximation to the yet non existing theory of quan- 
tum gravity. The range of validity goes down to the Planck-length, where quantum 
gravity effects are expected to become important. 

One of the main questions in this branch was concerned with the characterization 
of physical states on the various spacetime models. It is now accepted that all 
physical quantum states must have a singularity structure of Hadamard-form. For 
these states there is a well defined prescription to regularize the energy-momentum 
tensor and the semiclassical Einstein-equations 

Gab = 87r(T Qb ) 

which govern the backreaction of the quantized fields on the spacetime geometry 
make sense. Recently there was made progress in the formulation of the spectrum 
condition on a curved spacetime. There the Hadamard-states play an important 
role too, see ||. 

2. Dirac-field 

To set our notation we give in this section a short description of spinors in 
arbitrary curved spacetimes. Furthermore we recall the existence- and uniqueness- 
results of fundamental solutions of the Dirac-equation. At the end we mention 
Dimock's result that the Cauchy-Problem is well posed 0. 

By curved spacetime we mean a 4-dimensional C°°-manifold which is endowed 
with a metric, a covariant tensor-field of type (0,2), with Lorentzian signature 
(+,—,—,—). The metric describes locally the causal structure of our spacetime. 

l 
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We first introduce the complex Clifford-algebra of Dirac-matrices on Minkowski- 
spacetime and then we use the tetrad-formalism to lift them on our spacetime 
manifold. We denote the Minkowski-spacetime by (R 4 , r] a b — diag(l, —1, —1, —1)). 
Let C(R 4 ,?7) be the real, associative algebra with unit 1 which is generated by the 
elements {c(v),v e R 4 } and the relation 

{c{v),c{w)} = 2r](v,w)l (1) 

This is the Clifford-algebra associated to Minkowski-spacetime (R 4 ,?^). In the 
description of the Dirac-field on Minkowski-spacetime the so called Dirac-matrices 
play a prominent role. We obtain these matrices when we look at a faithful repre- 
sentation p of the complexified Clifford- algebra C(R 4 , 77) c on Mc(4). The resulting 
Dirac-matrices are denoted by j a :~ p(c(e a )), a = 0, ... 3. The vectors {e a } built a 
standard basis of R 4 . For two different representations p\ and p 2 of C(R 4 , ri)c we 
can always find a nonsingular matrix T which transforms the two corresponding 
sets of Dirac-matrices 7* and 7^ into each other, i.e. 

TjlT- 1 ^^ a = 0,... ,3 (2) 

In the following we assume a form of the Dirac-matrices with the additional property 

7o = 7o, l*j = -lj (3) 

The * means hermitian conjugation. Dirac-matrices with this property are said 
to be in standard representation. Next we consider the Lie-group Spin(l,3) of 
matrices S with the properties 

detS = 1 (4) 
SiaS- 1 = lb K\ (5) 

Together with the anticommutation relations of the Dirac-matrices one sees that the 
matrix A must be an element of the full Lorentz-group C. The mapping S — > A(S) 
defined by (5) is a 2-1 homomorphism of Spin(l, 3) onto C. In the following we re- 
strict ourselves to Spin (l, 3) the connected component of the identity of Spin(l, 3). 
The image under the above mapping is the group of proper orthochronous Lorentz 
transformations C\. 

Now let (M,g) the spacetime under consideration. We suppose (M,g) to be 
time- and space-orientable, which means that we can make a continuous distinction 
between future and past as well as left and right throughout M, and choose such 
orientations. Let TM denote the tangent bundle over M. To a given coordinate 
basis {e M , p = 0, . . . 3} at the point p e M we associate an orthonormal basis 
{e a , a = 0, ... 3} which is a linear combination of the coordinate basis vectors e M 
according to 

e a = e^ae^ (6) 
g(e a ,eb) = Vab (7) 
S^(p)eV"o > 0, e° > 0, p e M (8) 

The matrix (e M a ) shall be an element of GL(4, R) and the last condition is necessary 
to preserve the orientations. The collection of all this orthonormal bases constitute 
the so called orthonormal frame bundle F(M, g) associated to TM. Under a change 
of frame the tetrad e^ a transforms as 

e'\ - e\K\, A e 4 (9) 
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If we perform a local change of coordinates around peMwe obtain the following 
transformation law for the tetrad 



e 



(10) 



In this way F(M, g) admits the structure of a /^-principal fiberbundle. If we con- 
sider Spin (l,3) instead of C\ we obtain a Spin (l, 3)-principal fiberbundle which 
is called a spin structure S(M, g) for (M, g), if we have in addition a bundlehomo- 
morphism tp : S(M,g) — > F(M,g) which satisfies 

(poR s = i? A(s) o ip, S E Spin (l, 3) (11) 

The symbol R denotes the right action with A resp. S(A) on orthonormal basis 
vectors in F(M,g) resp. S(M,g). If a given spacetime possesses spin structures 
is determined by the topological properties of the underlying manifold. Here we 
consider only globally hyperbolic spacetimes for which the manifold is so well be- 
haved that the existence of spin structures is always guaranteed, but they are not 
uniquely determined by the topology. 

Now we are ready to define Dirac-spinors. On our spacetime (M,g) with spin 
structure (S(M,g),ip) we take a standard representation p of the Clifford-algebra 
C(R 4 ,?7)e together with the group Spino(l,3). The bundle of Dirac-spinors is 
now the associated vector bundle D p (M,g) to the Spino(l, 3)-principal fiberbundle 
S(M,g)), the representation p and the vectorspace C 4 : 

D p (M,g) := (S(M,g) x p C 4 )/Spin (l, 3) (12) 

In the following we abbreviate D p (M,g) with DM. A spinorfield is now a smooth 
section in DM. The collection of all smooth spinorficlds is denoted by C°°(DM). 
Those which have in addition compact support are denoted by C™(DM). With 
D*M we denote the dual spinorbundle. A v € D*M is for p E M an element of 
the dual space to the fiber over p. The smooth sections in D*M are called co- 
spinorficlds. The set of all smooth co-spinorfields are denoted by C°°(D*M). Let E 
be a submanifold of M. By DMs we denote the vectorbundle -K~j^ M (T 1 ) over E with 
projection Kdmyz '■= "\dm|tt2> m (E), where ttdm is the projection in DM. For a 
spinor u e C°°(DM) and a co-spinor v G C°°{D*M) we have a natural dual pairing 

v(u)\ p = v A u A , pe M (13) 

with respect to bases {E A , A = 0, ... 3} resp. {E A , A = 0, ... 3} in C°°(DM) resp. 
C°°(D*M). In the above expression we use the Einstein summation convention. 
We can now generate arbitrary spinor-tensor fields, i.e. smooth sections in the 
fiberwise tensor products of TM,T* M, DM, D* M . An element / in (<g) p TM) <g> 
(0 T*M) (g) (0 r DM) <g> ((g) s D*M) is given by specifying the family of C-valued 
functions / ai -°" ft h Al ^ Ar on M, such that 

•' 6l...0g B 1 ...B S ' 

/ = (f ai - ap bl ... bq Al - Ar Bi ^ B )e ai «...®e ap ®... 8 #«... 8 £ B ' (14) 

Of special interest is the spinor-tensor 7 <G T*M ® DM ® D*M which has the 
components 

la A B = {la) A B (15) 

with the Dirac-matrices j a , a — 0, . . . 3 in standard representation. For a vector 
field k E C°°{TM) we denote by jk the contraction of k with 7, i.e. jk A B := k a -y a A B . 
Next we introduce the notion of an adjoint spinor u + for a spinor u. This is 
the co-spinor with the following components: (u + ) B := u a jqab, where the bar 
means complex conjugation. We can now use the spinor-tensor 7 to lift the Levi- 
Civita-derivative of the metric g to mixed spinor-tensor fields in the following way. 
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The usual covariant derivative V : C°°{TM) -> C°°{T*M ® TM) is defined by 
(Vfc) b a := d a k b + T b ac k c , k € C°°(TAf). We extend this definition to spinor fields 
as follows: V : C°°(DM) -> C°°(T*M <g> DM) is defined through the specification 
of the components of V/, / G C°°(DM) with respect to bases (-E" 4 ) and (e„) as 
(V/) B := <9 a / B + o a B a! A i where the so called spin connection coefficients are 
given by 

*a B A == ~\Tca b l b B D l CD A (16) 

This defines indeed a covariant derivative on M. In demanding the Leibnitz-rule 
and commutativity with contractions, we extend this definition to all spinor-tensor 
fields. As an easy consequence of the above definitions we obtain the following 
lemma 

Lemma 1. The spinor-tensor 7 is covariant constant, i.e. V7 = 0. 

Proof. One can see this immediately by using the anticommutation relations of the 
Dirac-matrices ■ 

We can now write down the Dirac-equations for spinors and cospinors, which are 

{-ift + m)u = 0, u G C°°(DM) (17) 
(ift + m)v = 0, v€C°°(D*M) (18) 

In these equations m is a fixed real number. The Dirac-operator obeys Lichnerowicz' 
identity: 

(-« ft + m)(i ft + m)u = (n- -R + m 2 )u, u € C°°(DM) (19) 

where □ = rj ab W a Wb is the spinorial wave operator. The operator on the right 
side of (19) is known as the spinorial Klein-Gordon-operator. Next we consider the 
classical solutions to the Dirac-equation on a globally hyperbolic spacetime. We 
remind ourselves that a spacetime with given time- and space orientations is said 
to be globally hyperbolic, if it admits a Cauchy-surface. A Cauchy-surface is a 
smooth spacelike hypersurface ScM which is intersected by every endless causal 
curve exactly once. A globally hyperbolic spacetime has the structure RxE, i.e we 
can foliate M into smooth hypersurfaces. We now introduce distributions on the 
spaces C^(DM) and C°°(D*M). The natural pairing v(u) between v G C°°{D*M) 
and u £ C^°(DM) can be integrated over M in the obvious way: 

{v,u) M := / v(u){ P )dfi(p), u G C^{DM), v G C°°{D*M) (20) 

Here dfj,{p) is the volume element associated with the metric g. In addition to that 
we demand continuity of v, seen as a distribution over Cq°(DM), with respect to 
certain Sobolev- norms, sec JlTj for details. In this way we have an embedding of 
C ca (D*M) in Cq°(DM)', the space of cospinor- valued distributions over M. Anal- 
ogously {DM) is embedded in the space of spinor- valued distributions over M 
with compact support, C°°(D*M)' . The next Theorem gives us the existence of 
unique fundamental solutions to the Dirac operators 

Theorem 1 (Dimock, 1982). The operator (—ift + m) on C°°(DM) has unique 
retarded resp. advanced fundamental solutions : C§°(DM) — > C°°(DM) with 

{-ift + m)S ± = S ± (-ift + m) = id onC^(DM) (21) 
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and supp(S ± f) C J ± (swpp(/))J^] The operator {iJU + m) on C°°(D*A1) has unique 
retarded resp. advanced fundamental solutions S± : C^°(D*M) — » C°°(D* M) with 

{ip + m)S± = S±{ip + m) = id onC^(D*M) (22) 

and supp(S±f) C J ± (supp(f)). 

The difference S between retarded (S + ) and advanced (S~) fundamental solution 
is called Dirac-Propagator for spinors: S = S + — S~ . S$ = S+ — 5_ is called Dirac- 
Propagator for cospinors. 

The Cauchy problem for the Dirac-equation is now well posed: 

Theorem 2 (Dimock, 1982). Let uo G Cq°(DMs). Then there is a unique u G 
C°°(DM) with 

(—ift/ + m)u = and u|s = uq (23) 

In addition to that we have supp{u) C J + (supp(u )) U J~ (supp(u )). 

These two theorems lie at the heart of the quantization of the Dirac-field which 
is described in the next section. 



3. Quantization 

In this section the quantization of the Dirac-field on an arbitrary curved space- 
time is described. Since on a general curved spacetime we do not have time transla- 
tion symmetry, we have no natural criterion to single out a Hilbertspace of solutions, 
oscillating with positive frequency, which play the role of a one-particle space in 
Minkowski-spacetime. Therefore we adopt the algebraic approach here, in which 
only the algebraic structure of the theory is specified, i.e. we construct a C*-algebra 
and in a subsequent step certain subalgebras which are associated to open regions 



in spacetime |20| . These local algebras constitute the so called observable algebras 
in the sense of Haag and Kastler ||. A state is then specified as a positive nor- 
malized linear functional on the observables. The connection to the usual operator 
approach is given by the GNS-construction. 

Let S be the space of smooth complex-valued solutions to the Dirac-equation, 
which have compact support initial data: 

S = {ii£ C°°(DM)\(-ip + m)u = and u| s = w G C^(DM^)} (24) 

On S we introduce the following scalar product (.,.): S X S —► C, given by 

(ui,u 2 ) := / ui + (jhu 2 )(p)d^(p) (25) 

In this expression u is the restriction of the spinorfield u to S. n is the unit normal 
vectorfield on S and /i is the volume element induced by the three-dimensional 
Riemannian metric h a b = g a b + n a nb on E. This scalar product is independent of 
the choice of Cauchy-surface S. We complete the space S in the norm induced by 
this scalar product and obtain in this way a complex Hilbertspace H. Now we can 
associate to H an abstract C*-algebra .F[H] by the antilinear mapping B : H — > ^"[H] 
such that {B(h) : h G H} together with the unit 1 generate T[Y\] algebraically and 
fulfill the canonical anticommutation relations: 

{B(h),B(g)}=0 (26) 
{ B(h) , B(g)* } = (h, g) ■ 1, for all h,g G H (27) 



^^Here for N C M J ± (A r ) is the set of points in M, which are connected to TV through a future 
directed (+) resp. past directed (-) causal curve. 
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The algebra .F[H] has a unique C*-norm || . ||c*- Moreover it follows by the use of 
the CAR-relations that the mapping B is an isometry, i.e. we have 

IW)Hc = ||/||, / GH (28) 
By assigning to every open relative compact subset O C M the C*-subalgebra 
T{0) which is generated by elements {B(f) : f £ C^°(DM,0)}, where we denote 
by C^°(DM,0) spinors whose support lies in O, we obtain a net of field algebras 
O —> T{0). The even part hereof, i.e. elements of the form {B(h)*B(f) : h, f G 
Cq°(DM, O)} built the net of observable algebras O — > A(0) in the sense of Haag 
and Kastler f|. 

Now we come to the consideration of states on our CAR-algebra .F[H]. A state 
is a positive, normalized linear functional to on .F[H], by which we mean a linear 
map lo : J-[H] — > C with the properties 

uj(B(f)*B(f)) >0 (positivity) (29) 

= 1 (normalization) (30) 

As usual we focus attention on quasifree states, which have the form 

n 

u(B{h) . . . B(f 2n )) = sgn(a) JJ W (B(/„( 3) )B(/ ff(j+ „.))), n e N 

3=1 

(31) 

where the sum runs through all permutations a of {1, . . . , 2n} with cr(l) < cr(2) < 
. . . < a(n) and a(j) < cr(j+n), j = 1, . . . ,n. sgn(er) is +1(— 1), if er is an even (odd) 
permutation of {1, . . . , 2n}. In addition to that lo vanishes on all odd monomials 

B(/i)...%+i),«eN„. 

Let P be an orthogonal projection in H. A quasifree state up on J-[H] can then 
be specified in the following way. We choose an orthonormal basis {/i}igN for H 
consisting of eigenspinors for P. Then every element h € H has an expansion, 

h = EiEiNrC/i' h )fi- We denne an operator r on H by Th = ^i G N (/<> ^Z*- r is a 
well defined involution and it is antiunitary, i.e. 

(Tg,Th) = (h,g), Vg,h€H (32) 

Moreover T commutes with P. Then the spaces PH and (I — P)H together span 
the entire Hilbertspace H and are orthogonal w.r.t. (.,.). Therefore we can write 
every element ft € H in a unique way as 

h = h + + h~, h+ e PH, h~ e (l-P)H (33) 

Now we take the antisymmetric Fockspace J- a {\\) over H with annihilation op- 
erator a(f) : A" +1 H -> A"H defined by 

a(f)n - (34) 

and creation operator a(f)* : A™H — * A n+1 H given by 

o(/)*n = / (35) 

a(/)*(/i A . . . A /„) = / A f x A . . . A /„, V/,/i, . . . ,/ n G H (36) 

where fl is the vacuum vector in .F a (H). With A n H we denote the antisymmetric 
tensor product P A (<8>"H) defined by 

Pa(/i O . . . ® /„) = (n!)" 1 ^ xaM . <g) . . . ® / a(n) 

>?£S n 

where the sum runs through all permutations a in the permutation group S n and 
Xa is the sign of the permutation a. Furthermore we have used the notation 

/lA...A/ n = (n!)*P A (/i® ...®/„) 
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The operators a(f) and a(/)* are bounded operators on JF a (H) and they have 
anticommutation relations given by 

{a(f),a(gy} = (f,g)-l (37) 

This is a concrete C*-algebra, the so called Fock-representation of the CAR. We 
take the Dirac-field operator to be 

<M/) = a(f + ) + a(Tf-)*, f = /+ + /- where /+ € PH, f~ € (I — P)H (38) 

Then we have the anticommutation relation 

{*(/), *(<?)*} = (/,<?)•! 

This is known as the quasifree representation of the CAR-algebra .F[H] 
lop is the quasifree state determined by the two-point function 

Ar(/ ) 9) = (0,*(/)*( 5 )*0) 

where (. , .) is the inner product in P a (H). 

4. Hadamard states 

In this section we are looking at those states for which the two-point function 
has the singularity structure of Hadamard form ||. For these states it has been 
proven in the case of the Klein-Gordon field on a curved spacetime that the resulting 
operator description, which was sketched above, gives a physical theory in the sense 
that the principle of local definiteness is satiesfied |], |l6| . We define now Hadamard 
states for the Dirac-field which are believed to be equally important here as they 
are for the scalar case. For the introduction of these states we need the concept 
of bi-spinors and some technicalities concerning the causal structure of spacetime 
which we relegate in an appendix. 

Let (M,g) be a globally hyperbolic spacetime which Cauchy-surface E. The 
square root of the bi-spinor valued VanVleck-Morette determinant is defined by 

U A B ,=^/\q,q')J A B , (41) 

where A 1 ^ 2 (q,q') = g(q) 1 ^ 2 det(a(q,q') ;a b)g(q') 1 ^ 2 is the scalar VanVleck-Morette 
determinant and cr(q,q') is one half of the square of the geodesic distance s(q,q'). 
The bi-spinor of the parallel transport is defined by 

<r»J A B-# = (42) 

J A B >(q,q) = lV (43) 

Let V m A B i , to € No be a sequence of bi-spinors which are determined by the 
Hadamard recurrence relations Both U A b> and V m A B , are well defined on A, 
the set of points which can be joined by a unique geodesic (see the appendix for a 
precise definition of A). Furthermore we define for all n E N 

n 

V^ A B ,(q,q') := £ V m A B ,(q, q')(s(q, q')) m , (q,q') G A (44) 

m=0 

Finally we take a global time function T : M — > R and set 

rf («, q') := -2ie(T(q) - T{q')) - e 2 , q, q'eM,e>0 (45) 

When we look at the Dirac-operator we observe that it is not a hyperbolic oper- 
ator in the sense of Leray |l2j. We use the fact that, due to Lichnerowicz' iden- 
tity, the product of the Dirac-operator times it's complex conjugate is an hyper- 
bolic operator. In the next definition we will introduce the singularity structure of 
Hadamard form for the spinor Klein-Gordon operator which underlies the definition 
of Hadamard state for the Dirac-field. 



(39) 
. The state 

(40) 
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Definition 1. Let A be a sesquilinear form on C^(DM). A has Hadamard form 
(for the spinor Klein-Gordon operator), if the following data exist: 

1. a causal normal neighborhood ./V of a Cauchy-surface £ 

2. a .ZV-regularising function \ 

3. a smooth future directed global time function T on M 

4. a sequence H (n *> G C n (DN K D*N), n G N, so that for all n G N and all 
f,h€C$°(DM,N) 

A(f,h) = lhn| (A^' n ) A B ,(q,q')f^(q)h B ' (q')dfj,(q)dfi(q'), where 
(Af'")\, (?,</) :-x(<Z,9')(Gr'") A s ,(g,g') + ^ (ri) V(g,<z') with 

(G^Y* (?,</) := ^(-^^^ + V^ A B ,( q , q ')H S ( q , q ') 
Aiz l \s{q,q')+ri(q,q') 

+ rj( q , q ')) 

for (q, q') G A n (AT x JV); In is the main branch of the logarithm 

We now view the two-point function u)(B(fi)* B(f2)) as a hermitian form /i, /a i— > 
(/ijQ/2) over H. This form is a solution of the Dirac equation in both arguments, 
i.e. we have 

(H A 7 + m)fi,Qf 2 ) = (h,QHN + m)f 2 ) = (46) 
With the next definition we arrive at the Hadamard form of a state of the Dirac- field 

Definition 2. A quasifree state u> on .F[H] is called a Hadamard state of the Dirac- 
field, if there is a solution A of the spinorial Klein-Gordon equation of Hadamard 
form, so that the two-point function of u> has the following form 

w(B(A)*B(/ 2 )) = A((t/7 + m)/!,/ 2 ), /i,/ 2 G H (47) 

A is a solution of the Dirac equation in the first argument, since it fulfills Lich- 
nerowicz' identity and in the second argument because of hermiticity. The singular 
part of A is determined by geometric quantities like metric and spinor connection 
only. Therefore a state is characterized, if we fix the smooth part of A. 



Since the work of Radzikowski 1 15 there is another possibility to characterize the 
singularity structure of a Hadamard state. This method uses Hormanders wavefront 
sets to describe the singular behavior of the numerical distributions 

/1, ■ ■ ■ ,/n H A£(/i, ... , /„) := (fi, *(/!)... *(/„)«> (48) 

We will here only give a brief introduction to this powerful theory; for more informa- 
tion see p], H, |]. The idea in using wavefront sets is to localize the distributions 



around the singular support and then to analyze the directions in Fourierspace 
which causes these singularities. The advantage of this approach is that only local 
concepts are used. In the following we will give the definition of the wavefront set 
of a distribution and recall the main result of Radzikowski's thesis. 

Definition 3. Let (M, g) an n-dimensional manifold with associated cotangent 
bundle T*M. Let u G D'{M). A point (x ,k ) € T*M \ {0} is called a regular 
directed point of u G D'(M) if and only if for all Ao G R n and for every function 
4> G C°°(M x R™, R) with d x (j>(xQ, Ao) = ko there exists a neighborhood U of xq in 
M and a neighborhood A of A in R", so that for all p G C°°(U) and all TV > 
uniformly in A G A: 

\{u,pexp-* T ^-' x >)\=o(r- N ) for r -> 00 (49) 



ADIABATIC VACUUMSTATES OF THE DIRAC-FIELD ON A CURVED SPACETIME 9 



The wavefront set WF(u) is then the complement in T*M\ {0} of the set of regular 
directed points of u. 

The wavefront set of a distribution u has some remarkable properties: 

1. The projection of WF(u) on the first factor gives singsupp(u) 

2. WF(m) is a closed subset of T*M \ {0}, since the complement of WF(u) are 
open sets. 

3. For all testfunctions <fi with compact support it holds 

WF(0u) c WF(tt) 

4. For a differential operator D with C°° coefficients we have 

WF(Dh) c WF(u) 

The last relation means that the application of a linear differential operator D on 
u does not enlarge the wavefront set of u. 

For the wavefront set of a tensor product u ® v of two distributions u and v on 
M the following formula holds (see Hormander ||): 

WF(u ® v) C WF(ti) x WF(«) 

U ((supp(u) x {0}) x WF(u)) (50) 
U (WF(tt) x (supp(v) x {0})) 

With the help of wavefront sets we are now able to characterize the Hadamard 
states of the Dirac field. For this purpose we first look at the wavefront set of 
the two-point function of the Klein-Gordon field. As was explained previously we 
obtain the two-point function of the Dirac field through the application of the 
adjoint Dirac-operator on the scalar two-point function. Due to the fact that this 
procedure does not enlarge the wavefront set (since the Dirac-operator is a linear 
differential operator) it is sufficient to determine the wavefront set of the two-point 
function of the Klein-Gordon field. A completely local description of Hadamard 
states of the Klein-Gordon field was given in |l5| : 

Theorem 3 (Radzikowski, 1992). A quasifree state co of the Klein- Gordon field 
on a globally hyperbolic spacetime (M,g) is a Hadamard state, if and only if it's 
two-point function Aj K q has the following wavefront set: 

WF{Kl KG ) = {(n.fi; z 2j -6) G T*M 2 \ {0}; (n.fc) ~ 0r 2) £ 2 ),£° > 0} (51) 

Here (xi,£i) ~ (^2,6) means that x\ and X2 are connected by a lightlike geodesic 
7 in the way that is tangential to 7 in x% and £2 is £1 parallely transported 
to x 2 . On the diagonal (xi,£i) ~ (2^2, £2) means that £1 = £ 2 and (£i) 2 = 0. We 
obtain the following result as a corollary to the above theorem: 

Corollary 1. A state u> on ^"[H] is a Hadamard state for the Dirac-field, if and 
only if it 's two-point function A£ has the following wavefront set 

WF(A%) = {(n,ei;a; 2 ,-6) G T*M 2 \ {0}; ( Xl ,^) ~ (x 2 ,b),$ > 0} (52) 

Proof. This is an immediate consequence of property 4. above. I 

One finds the wavefront set of an arbitrary m-point function with the help of (50). 
When we denote the kernel of A!^ by A^(xi, ... , x m ) we obtain 

n 

WF(A£ (an,... ,x m ))c(J0WF(A^' ) ), n = m/2 (53) 

a j = l 
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where is the two-point function in the variables aWj) > x a (j+n) viewed as a 

distribution over M m , with wavefront set 

WF(A2 (j) ) = {(xi,0;... ;x a(j) ,k aU) ;... ; x a(j+n) , k a(j+n) ; . . . ;x m ,0)\ 

,.s (54) 

(^OJ^aOJ^aO+nJ^^O+n)) € WF(A2 Uj )} 

That this expression is the wavefront set of A^ is clear, since A^ is a sum over 
tensor products of two-point functions. 



5. DlRAC-EQUATION ON ROBERTSON- WALKER SPACES 

In this section we consider homogeneous and isotropic spacetimes. It is known 
that the requirement of homogeneity and isotropy leads to the Robertson- Walker 
spacetimes of the form M — R x £ K , where S K is a homogeneous Riemannian 
manifold. The parameter k can take the values +1,0 or — 1. Correspondingly 
S K is a space of constant positive, vanishing or constant negative curvature. We 
focus attention to the Dirac-equation for which we construct explicit solutions by 
the method of separation of variables due to Villalba and Percoco OL We will 
use these solutions in the next section to define adiabatic vacuum states of the 
Dirac-field in this spacetime. 

The metric in chronometrical coordinates r, \ , 6, if takes the form 

ds 2 = c a{T \dT 2 - d X 2 - £, 2 {x)dn 2 ] (55) 



where dfl 2 is given by 

dn 2 = d0 2 + sin 2 9dip 2 (56) 



(6 £ [0,7r], ip £ [0, 2n], x e [0, oo) for k = 0, -1 and x € [0, n) for k = +1) and 
a(r) is a real- valued smooth function. The function £(x) is given by 

fsinx falls k = +1 
X falls k = (57) 
sinhx falls K = — 1 

To calculate the spinor connection we use a vierbein basis e a ^ defined by 

efelnab = <V (58) 
The formula for the Christoffel symbols in this basis is given by 

r\ c = e *e£ (a M e v c + e c A r% A ) (59) 



In this expression T v are the Christoffel-symbols in a coordinate basis. Greek 
letters will always denote indices with respect to a coordinate basis, whereas Latin 
indices correspond to a vierbein basis. We choose the vierbein basis in the following 
way: 



p a = e «/2( S aO Sti0 + galg^ + fa^S^ + fa) jm^) (60) 

e v b = e- a / 2 (6 u0 S b0 + 6 ul 6 bl + C l {xW 2 h-2 + CHx) sin^ 1 0<T 3 <5 63 ) (61) 
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With this input one finds after a straightforward computation the following expres- 
sion for the Christoffel symbols in vierbein basis: 



1 



a ' 2 \a5 w 



+ 5 bl 
+ 5 b2 

+ S b3 



<5al<5 c l + 5a25 c 2 + 5 a3 S c3 

S al 5 c0 a - 25 a2 5 c2 C\x)t'(x) ~ 26^5^ (x)C \x) 
S a2 S cQ a + 25 a2 5 cl C 1 (x)Z'(x) ~ ^S^S^Hx) cot# 
S a3 S c0 a + 25 a3 6 cl C 1 (x)Z'(x) + 26 a3 6 c2 C 1 (x) cotO 



(62) 



The spinor connection is given by 



°a A ■= ca 7& D l A 



(63) 



The 7 a are the components of the spinor-tensor 7 w.r.t. the used moving frame. 
The explicit expression for the spinor connection coefficients can now be given: 



1 



-a/2 



^al707 1 + <>a2707 2 + ^a3707 3 

2s a2lll 2 r 1 (x)C(x) ~ 2£ Q37l 7 3 r 1 (x)£'(x) 

2^3727 3 r 1 (x) cot/ 



(64) 



Here the anticommutation relations of the Dirac-matrices were used. 
The Dirac-equation reads 

H 7 a V Q + m)^{x) = (-i'fe£V ll + m)^{x) = 



(65) 



Here the covariant derivative is lifted to the spinor bundle DM, as was explained 
previously: 



V Q = d a + a a 



(66) 



Inserting the expressions for the spinor connection coefficients and the vierbeins we 
obtain the following explicit form of the Dirac-equation: 



-ic 



-a/2 



7°(9o-^)+7 1 (ai-r 1 we , (x))+ 



l 2 r 1 (x)(d 2 --cote) + 7 3 r 1 (x)sin- 1 05 3 ] +mU(x) =0 



(67) 



In order to arrive at a more simpler equation we make the following transition to 
the new spinor <3> defined by: 

$(a;) ^C 1 (x.)sm- 1/2 9e~ 3a/4 ^(x), x ^ 0, 8 ^ 0, tt, 2tt (68) 

The Dirac-equation for $ now reads: 



7°do+7 1 di+7 2 r 1 (x)d2+7 3 r 1 (x)sin- 1 0d 3 + m = (69) 



-a/2 



We now want to solve this equation by the method of separation of variables. 
This was already done in the paper by Villalba and Percoco jl^] , so we shall outline 
here the main steps of this calculation only. All of the missing details are completely 
straightforward. 
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First of all we decompose the Dirac-operator into two commuting differential 
operators which can be dealt with separately. We use the Minkowski-space Dirac- 
matrices in the following representation: 

Mo -u) «)■ ™ 

The Pauli-spinmatrices a* are given by 



1 OJ \i J \0 -1, 

The Dirac-equation can be written in the form: 

[D 1 (e^) + Dl(T, X )]$(x)=0 (70) 

Here we have defined $ = 7 1 7°$ and the operators D 1 (9, tp) and D 2 (t, \) are given 
by: 

D\9, tp) = [ 7 2 a 2 + 7 3 sin" 1 9d 3 ] 7 1 7 ° (71) 

Dl{r, x) = ^(X) [l°do + j'di + me"/ 2 ] 7 1 7 ° (72) 

We observe that the operator D 1 (9,p) is sclfadjoint with respect to the scalar 
product: 

(*,$):=/ Ml d^ A {6,tp)+$ A {0,tp) (73) 
Jo Jo 

Now we have the following lemma: 

Lemma 2. D 1 and D 2 are commuting operators on C yo (DM), i.e. we have 

[D\6, i p),D 2 K (r,x)]=0 (74) 
Proof. We evaluate the product D 1 D 2 : 

B 1 Bl = [7V7 3 2 +7 3 7 1 7°sin- 1 ^ 3 ]£(x)[7 do + 7^1 +™e a/2 ]7 1 7° 
= 7 1 7° b 2 d 2 + 7 3 sin" 1 9d 3 ] £( X ) [7% + 7% + ime a/2 ] 7V 
= - 7 1 7°£(x) [l°do + 7^1 - ime a ' 2 } [ 7 2 5 2 + 7 3 sin" 1 flflbfrV 
= £(x) [j°do + 7^1 + ime"/ 2 ] 7 1 7 ° [ 7 2 5 2 + 7 3 sin" 1 0%] 7 1 7° 
= D?D 1 



Because of this lemma it should be possible to find a common set of eigenspinors for 
these two operators. Therefore we make the following ansatz to obtain a solution 
to equation (70): 

D\9, <p)*(x) = A$(x) = -D 2 K (r, xMx), A e C (75) 

First we look at the angle dependent part D 1 (9, tp): 

D l (9, tp)$ = A$ |7 2 a 2 +7 3 sin- 1 6»a3|7 1 7 $-A$ = (76) 

A solution to this equation can be obtained in terms of orthogonal polynomials, i.e. 
we have the result 

Proposition 1. 1. The eigenf unctions of the operator D 1 have the following 
form: 

$(e,tp) = U- 1 Z(6, V >) (77) 
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G n i(cosfl) = - 7 ^=P,° Jl,/Ja) (cosfl) (81) 

F nl (c OS e) = (-iy^=p^\ cos e) (82) 



where S = is a spinor with components 

— > = (^>) ™ 

wit/i /gN, n = |(2m + 1), m € N and scalar valued functions 

<#,(0,¥>) = sm n 6 cos(^/2)G„ ( (cos 0)W n (<p) (79) 
C£,(0,p) = sm"0sin(0/2)i^cos0)W„( V ) (80) 

G n i and F n i are the Jacobi-polynomials P^ a,b ^ in the variable cos0 of order I, 
i.e. 

1 

1 

The parameters [3\ and (32 are given by 

and N n i is a normalization constant. The scalar valued function W n is given 
by 

W n (<p) = exp(tnv») (83) 
And the unitary matrix U is given by 

U= i (7 V + 7V + 7V + 1 4 ) (84) 
2. The eigenvalues A are given by 

A = / + n+ieN\{0} (85) 
The Jacobi-polynomials are defined in the proof. 

Proof. Eq.(76) is a system of partial differential equations in two variables and the 
aim of the following calculation is to separate these two variables. The matrix U 
transforms the Dirac-matrices in the following way: 

U^U' 1 = 7 , U^U' 1 = 7 3 , Uj 2 !/- 1 = 7 1 , U^U' 1 = 7 2 

The application of U to equation (76) then gives 

{7 1 <92+7 2 sur 1 <9<9 3 }7 3 7 S-AS = (86) 

Inserting the above form of the Dirac-matrices we obtain 

+ "O*- 1 "*}©-^- ,87 ' 

The variable ip enters this system of equations in the derivative operators only. 
Therefore we make the following ansatz for the spinor S: 

Si=ai(r)6i( X )^^)wM, S 2 = a 3 (r)Mx) (*^) W(<p) 

This is the most general ansatz for the spinor S, for which one can separate the 
variables 9 and ip in eq. (87) . This form is also compatible with the structure of the 
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operator D^(t, X ) as we will see later. When we insert this into equation (87), we 
see that we have 

Xl = X3 X2 = -X4 

so we obtain only two coupled linear independent partial differential equations in 
the variables and ip. Now we can separate this system in the following way: 

sinflxJ^iAxi - d eX 2) = n = -iW^d^W (88) 

sm0 X ^(tXx2 + dexi) = n = -iW-'O^W (89) 

The solution for W reads: 

W n (tp) = exp(inp) (90) 

We demand the spinor 'J to change the sign when rotated about an angle 2tt in the 
direction tp, so n varies in the range n = ±5, ±|, ±|, .... 

Next we have to solve the remaining system of equations in the variable 0: 



(do + -^jX2 - Axi 

(de- ^-^)xi + Ax2 
V smu/ 



(91) 
(92) 



with A := iX. We make for xi an d X2 the ansatz 

X i„(0)=sin"0cos(0/2) 5 (0) (93) 

X2 „(0)=sin"0sin(0/2)/(0) (94) 

where / and g are unknown scalar functions, still to be determined. Inserting this 
into equation (92) we obtain after a little manipulation 



cos(0/2) 



ng(6) (cos 6 - l) - sin 2 0.g(0) 
sin(0/2)[sin0(A/(0)-i ff (0)) 



where 



With the identity 



we finally obtain 



g'(0) 



dg d cos 9 
d cos 9 d9 



-g(9) sin 9 



cot(0/2) = 



1 + cos ( 



sm( 



-ng{9) - (1 + cos0).g(0) + A/(0) - -g{9) = 



With equation (91) we proceed in the same way and find 

nf(9) - (1 - cos0)/(0) + l -f{9) - \g{9) = 

With x := cos0 these two equations become 

. dG (I 
{1 + x) ^ + [-2 +n 
* dF ( 1 

The functions F and G are defined by 

F(x) := f(9) 
G(x) := 5 (0) 



G = Ai^ 
F = AG 



(95) 

(96) 

(97) 

(98) 

(99) 

(100) 
(101) 

(102) 
(103) 
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When we now insert eq.(lOO) into eq.(101) we obtain a differential equation which 
is solved by the Jacobi-polynomials jy : 

(1 - x 2 )G + (fa -fa- (fa + fa + 2)x)G + 1(1 + fa + fa + 1)G = (104) 

Here we have made the definitions: 

ft -=n-\ 
fa:=n+- 

Z:=|A|-|n|-i 

The Jacobi-polynomials of order I are defined by 

i 

P[ 01 ' 02) (x) = 2- 1 Y, c m (fa,fa)g m (x\ I £ N where 

m=0 

I + fa\ (I + /V 



C m (fa,fa) = 

V lit 

g m (x) = {x-l) l - m {x + l) m 
The parameters fa , fa are restricted to the values 

fa > -1, fa > -1 

This restrict the values of the parameter n further to n = h (2m + 1), m G N The 
solution for G now reads 

G nl (x) = -1-P^\ X ) (105) 



where N n i is a normalization constant. The eigenvalues |A|, which are related to 
the order of the Jacobi-polynomials I and to n via I = |A| — n — i must be integer, 
since I is an integer. With G we easily obtain the solution for F: 

F nl (x) = ^=(-l) l P^\x) (106) 



If we choose N n i to be 

_ 2 2 "+ 3 7rr(/ + n + l/2)r(/ + n + 3/2) 
nl ~ (2l + 2n + l)T(l + 2n + l)U { ' 

then the eigenfunctions $(#, ip) of D 1 are orthonormal, i.e. we have the relation: 

ds f d^+(e,cp) A ^ n , v (e,ip) A = 5 nn ,5 w (los) 

Jo 

We now consider the operator D 2 (t, %). Equation (75) reads 

D 2 k (t,x)$ + \$ = 
7°<9 + 7 1 3i + ime a/2 \ 7 1 7 °l> + A^ x $ = 
After the application of U, this equation becomes 

-7 3 <9 - 7°<9i + 7 3 7°jme Q / 2 + A^ 1 ] S = (109) 
Inserting the Dirac-matrices we obtain the following system of equations: 
-a 3 d a E 2 ~ d{Zi - a 3 ime^ 2 E 2 + AT^i = 

/, - i (no) 



o 
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(111) 

(112) 



As in the angle dependent case we are looking for a separation of the variables r 
and if. Therefore we enlarge our ansatz for the spinor 2 in the following way: 

'ai(-r)&i(x)<^(^r 

as^baixjC^ip) 

-<H(r)b A (x)C*,V,<Ph 
where ai,bj, i,j = l,...,A are unknown scalar valued functions, to be determined 
in the following. As was mentioned earlier in the proof of proposition 1, to be 
compatible with both operators D 1 and D 2 we must have 

a\ = a 2 , b\ = b 2 
a a = a 4 , b 3 = 6 4 

This is indeed the case, as can be seen when eqs.(lll), (112) are inserted in (110). 
Therefore we have to solve only the reduced system 

a^daax + b^d^ - ime a/2 a 1 a^ 1 + A^" 1 ^ 1 = (113) 

-a^doas - b^dxh - ime a/2 a 3 a^ 1 + X^hb^ 1 = (114) 

Introducing the two constants z±, z 2 € C we obtain for the variable r: 



ime o/2 )ai(r) = Zia 3 (r) (115) 



(— 

And for the variable x : 

d 



ime a / 2 y 3 (T) = z 2 ax{T) (116) 



(^ + k- 1 (x))b 3 (x)=-zMx) (H7) 

(-^-k-Hx))bi(x) = -Z2bs(x) (US) 

The equations in the radial coordinate x can be solved by hypergeometric func- 
tions. 

Proposition 2. For n — the solutions to the eqn.(117) and (118) can be ex- 
pressed through the functions 

h( X ) = B lw (y) = B k {d x y\ k {y) + d 2 yy k {y)) (119) 

h(x)=B 3w (y)=B k (d 1 yj k - 1 (y) + d 2 yn-i(y)) (120) 

where d\,d 2 € R, k € Z, y — ±^fwx, with w € R+ 

Proof. First of all we take k = in (117) and (118) and combine them to the single 
equation of second order 

b'l - h {\ 2 x~ 2 ~ \X- 2 + z)=0 (121) 

where z = z\z 2 and a prime denotes a derivative with respect to x- F° r z S R := 
{z £ R|z < O}^ we set — z =: w > and define a new variable y by y := zt^/wx- 
The derivative w.r.t y is denoted by a dot. Then equation (121) becomes 

y 2 B 1 (v) + Bi(y)[X(l — A) + y 2 ] = (122) 



2 On physical grounds there are only negative values for 2 of interest, since it follows from the 
treatment of the equations (115) and (116) that —z is the absolute value squared of the wavevector 
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where the function B\ is defined by B\(y) :— b\(x)- If we finally set A = — k = 
—k(l,n) we obtain the so called Riccati-Bessel differential equation [0: 

y 2 ^ (y) + B 1 (y) [y 2 - k(k + 1)] = (123) 

The solutions to this equation are the spherical Bessel-functions of the first and sec- 
ond kind jfc(y) and yfc(y) respectively. The general solution is a linear combination 
of these two: 

B lw (y) = B k (d 1 yi k {y) + d 2 yyk(y)) , k € Z, d x , d 2 € R (124) 

where B k is a normalization constant. The solution for b% w is, expressed through 
the function B 3w (y) := 63™ (x), 

s 3 ™(y) = B ^y~ 1/2 ( d i 4-i/2 (v) +d 2 y k -i/2(y)) (125) 

where the Jfc+i/2(2/) and Y fc+1 / 2 (2/) are the usual Bessel-functions of rational order 
of first resp. second kind. They are related to the spherical ones by 

h(y) = J^y- 1/2 h+iMy), y*(y) = jr 1/2 Y fe+1/2 ( y ) (126) 



We summarize the above calculation in the following two statements 

Proposition 3. The solution to the Dirac- equation eq.(67) has the following spinor 
structure: 

V wn i (x) = - e 3Q / 4 X sin 1/2 9ME wnl (x) (127) 

where M is the matrix M = 7°7 1 [/~ 1 and B, wn i(x) is given by the equations (111) 
and (112). 

Proposition 4. The spatial part of the solution of the Dirac- equation (67) in the 
coordinates x, 0,tp on a fiat Robertson- Walker spacetime has the form: 



1. 



= j dw^(^a n i(w)^ wn i(x) +a* nl (w)^ wn i(x) S j (128) 
with basis f w „i given by 



^wni(x) = -xsin 1/2 9ME wnl (x) 

where 
with 



^wnl (^) 



^lwnl(x) = b lw (x) 



The coefficient function a n i(w) in eq.(128) must be chosen with compact support 
to make the integral convergent. The solution to the equations (117) and (118) in 
the cases n = ±1 can be given equally in terms of more advanced hypergeometric 
functions but there is only little to learn of these formulas, so we skip them here. 

The treatment of the equations (115) and (116) requires the specification of the 
function a, which corresponds to choosing a certain universe model. 
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6. Adiabatic vacuum states 

We come now to the discussion of adiabatic vacuum states on .F[H]. Associated 
with each homogeneous space E K there is a group of isometries G K . This group acts 
through g(t,x) — (t,gx), g € G K on M. In the case of constant positive curvature 
(k = +1) this is the rotation group SO (4). If we take flat Cauchy-surfaces we have 
G K = E(3). In the case of constant negative curvature (n = —1) the symmetries 
are the proper orthochronous Lorentz-transformations C\_. See ref. p3| for more 
details. 

We now want to characterize certain states on the CAR- Algebra .F[H] which are 
invariant under the action of the transformation group G K . The unitary represen- 
tation U of G K on C$°(DM), defined by 

(U(g)f)(x) := fig-'x), g E G K (129) 

commutes with the Dirac-propagator for spinors S : C£j°(DM) — > C°°(DM). In 
addition to that it holds that the operators U (g) commute with the antiunitary 
involution T on H. In this way they form a group of Bogoliubov-transformations 
on H. Induced by these transformations there are automorphisms a g oi J-[H], i.e. 

a g (B(f)) = B(U(g)f), geG K (130) 

These so called Bogoliubov-automorphisms are now used to single out homogeneous 
and isotropic states of the Dirac-field. A state ui on .F[H] is said to be homogeneous 
and isotropic, if and only if it is invariant under all Bogoliubov-automorphisms a g , 
i.e. 

Luoa g = u, V.g S G K (131) 

The restriction to quasifree states makes it possible to formulate an equivalent 
condition on the two-point function 

AS(/i,/j)-wWiWj)) (132) 

Definition 4. A quasifree state cu on .F[H] is called homogeneous and isotropic, if 
and only if it's two-point function satisfies 

^(U(g)fi,U(g)f 2 )=A%(fi,f 2 ), h,f 2 G H, g e G K (133) 

We now introduce adiabatic vacuum states for the Dirac-field on a Robertson- 
Walker spacetime. For the Klein-Gordon-field these states were invented by Liiders 
and Roberts |l3| , motivated by Parkers requirement that the particle production by 
the expanding universe should be minimal [ jj"4| . It has been shown that adiabatic 
vacuum states for the Klein-Gordon-field are the physically correct states, in the 
sense that they have the singularity structure of Hadamard-form ||. This fact 
motivates the same procedure in the case of the Dirac-field. First of all we write 
the field operator according to proposition 4 as: 




(134) 



where a n i(w), a n i(w) + are annihilation and creation operators on the antisymmetric 
Fock-space over H. They have anticommutation relations 

{a n i(w), a n >t>(w')*} = S(w -w')5 m >8u> (135) 

We now consider the equations (115) and (116) for the time dependent part. Com- 
bining them we get 

" *f me ' 2 + m2eQ ) ai -( T ) = ~ wa ^(r) (136) 
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In analogy to the Klein-Gordon case we make a WKB-type ansatz for a solution to 
equation (136) 



a lw {r) = e~ 3Q / 4 (2Cl w ( T )y 1/2 exp \i f dr'Sl^r') 

. Jo 



(137) 



where Q w is a complex valued smooth function with positive imaginary part. If we 
insert this in (136) we get an equation for Q w : 



9. 2 3.fUr) 3i l/fi™(T)\2 3.. 



(138) 



dme Q / 2 + m 2 e a + w 



2CI w (t) 2n w {r) 2 2 

In a slowly varying universe we expect the derivative terms to be small, so we try 
an iterative solution to this equation: 

(0(°)(r)) 2 ^ + raV 
(^ +1) (r)) 2 = - + mV - l -ame^ + ^ - + 3 -a^l 

' il -n(r) ( \ , 1(^W\ 2 l^ r) (r) , 1 f4 r) (r) 



The only freedom we have in the choice of a homogeneous and isotropic quasifree 
state is to specify initial data for the functions a\ w and a 3w . Therefore we make 
the following definition: 

Definition 5. An adiabatic vacuum state of the Dirac-field of order r on a flat 
Robertson- Walker spacetime is specified by the following initial data for the function 
a\ w : 

ai w (r) = W^(t) (139) 
ai w (r) = W£\t), where (140) 



Wi r \r) :=c- 3 «/ 4 (2^(r))- 1/2 exp L frf/fiW^ 

. Jo 



') 



(141) 



The function as w is then fixed by eq.(115). Now one calculates the two-point 
distribution to be 

A 2 (0:1,0:2) = (f2,*(a;i)*(a;2)fi) 

dw^2^ wn i(xi)y wn i(x 2 ) 



I- 



where the anticommutation relations eq.(135) were used. To decide weather the 
quasifree state associated with this two-point distribution is an Hadamard-state we 
must look at the short distance behaviour of the expression (142). If we use the 
Taylor expansions of the Bessel-functions we arrive at a form which agrees with the 
one in definition 1. So we conjecture that we obtain in this way a Hadamard-state. 

7. Appendix 

In this appendix we define bi-spinors in curved spacetime and introduce the rele- 
vant notions of the causal structure of this spacetime for the definition of Hadamard 
states for the Dirac-field. In the following our spacetime is a smooth C°°-manifold 
on which there is defined a Lorentz metric g. In addition to that we assume that 
(M, g) is globally hyperbolic. 
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Let A be a open subset of M. Further let DNMD* N be the outer tensor product 
of the bundles DN and D*N. This is the bundle over the product manifold A x A. 
the fibres at (p,p') G TV x TV are (DN) p ® (D*N) p /. ir is the product-projection, 
i.e. 7r((p, /), (p', /')) = (p,p')- Smooth sections in DN^D*N are called bi-spinors. 
With the choice of vierbeins (25a) a t P and (2?' s ) at p' , we construct a new vierbein 
(q,q') i ► (2?H2?')A B '(g,9') := E A {q) ® E' B ' {q'), A, B' = 0, . . . , 3 in 7JA HI D* A 
at (p,p')- Smooth sections v € C°°(DN S3 Z?*A) are given by 

ufe g') = « V (g, & ^')a S ' (g, g') (143) 

with C°°-functions u^s/. 

For the singularity structure of the two-point function to be well defined we 
need the notion of a causal normal neighborhood of a Cauchy-surface E, which we 
explain now. An open neighborhood A of a Cauchy-surface E in (M, q) is called a 
causal normal neighborhood for E, if for every choice of p, q G TV with p G </ + (g) 
there is a convex normal neighborhood pll 7?. C M, so that J~(p) H J + (q) C 7?.. 
It was shown by Kay and Wald jl0| that to every Cauchy-surface there exists a 
causal normal neighborhood. Now let A be the set of points (p, q) G M x M, such 
that p and g can be joined by a causal curve with J + (p) H J~(q), if g G «2 + (p) 
or J~(p) H J + (q), if p G J + (q) contained in a convex normal neighborhood. The 
square of the geodesic distance s(p, g) of p and g is then well defined and smooth 
in X C M x M. 

Let TV be a causal normal neighborhood of a Cauchy-surface S. A smooth 
function x on TV x A is called a A-regularizing function, if it has the following 
properties: There is an open neighborhood Y C N x A of the set of causally 
connected points in A, such that Y C X and x = 1 on 7 and \ = outside of A. 
Then x is well defined, since Y and the complement of A in A x A are disjoined 
closed subsets of A x A. It can be shown that such a A-regularizing function 
always exists (l5). 
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